DCPT-02/43 

Bound monopoles in Brans-Dicke theory 



Betti Hartmann* 

Department of Mathematical Sciences, University of Durham, Durham DH1 3LE, U.K. 

(February 1, 2008) 



Abstract 

We consider axially symmetric SU(2) Yang-Mills-Higgs (YMH) multi- 
monopoles in Brans-Dicke theory for winding number n > 1. In analogy 
to the spherically symmetric n = 1 solutions, we find that the axially sym- 
metric solutions exist for higher values of the gravitational coupling than in 
the pure Einstein gravity case. For large values of the gravitational coupling, 
the solutions collapse to form a black hole which outside the horizon can be de- 
scribed by an extremal Reissner-Nordstrom solution. Similarly as in the pure 
Einstein gravity case, like-charged monopoles reside in an attractive phase in 
a limited domain of parameter space. However, we find that the strength of 
attraction is decreasing for decreasing Brans-Dicke parameter to. 
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I. INTRODUCTION 



Topological defects |T| are believed to have formed during phase transitions in the early 
universe and are related to spontaneously broken symmetries. Depending on the topology 
of the vacuum manifold, <i-dimensional objects form. The d = 0, 1, 2 defect is the magnetic 
monopole, the string and the domain wall, respectively. Two basic types of topological 
defects have been considered in the literature : a) defects which are described by theories 
with a broken global symmetry and b) defects which are described by theories with a broken 
local, i.e. gauge symmetry such as the Nielsen- Olesen vortex || and the 'tHooft-Polyakov 
monopole ||. While the latter exhibit particle-like properties such as finite energy and a well 
defined core, the former have divergent energies. In turn, due to their long-range behaviour, 
global defects have much stronger gravitational effects than their local counterparts 

MM. 



The local magnetic monopole results from a spontaneous symmetry breaking of an SU(2) 
gauge symmetry down to a U(l) gauge symmetry. The mapping of spatial infinity to the 
vacuum manifold is characterized by the winding number n. The spherically symmetric 
n — 1 solution is the famous 'tHooft-Polyakov monopole [BJ. Since this was proved to be the 
unique spherically symmetric solution ||, the construction of multimonopoles longed for an 
axially symmetric Ansatz which was introduced in 



In flat space like-charged monopoles are non-interacting |TTJ in the limit of vanishing 
Higgs boson mass. In this so-called Bogomol'nyi-Prasad-Sommerfield (BPS) [|l|,fl limit, 
the attraction of the long-range Higgs field exactly compensates the repulsion caused by the 
long-range U(l) field. These configurations saturate a lower energy bound, the so-called 
Bogomol'nyi bound such that their energy is proportional to the winding number n. For 
finite Higgs mass, the monopoles are repelling ][H| since now the Higgs field is exponentially 
decaying. 



Gravitating monopoles have been studied in the framework of General Relativity fL4,15 



Interestingly, it was found that gravitating monopoles can - for specific choices of the gravi- 
tational constant and the Higgs self-coupling constant - reside in an attractive phase in the 
sense that their energy per winding number is decreasing with increasing n |L5| . 



Among alternative gravity theories, the scalar-tensor theory introduced by Brans and 
Dicke [16] is one of the most popular. Since Mach's principle ]L7|] leads to problems within 
the framework of General Relativity, Brans and Dicke introduced a scalar field which plays 
an analog role than the inverse of Newton's constant and is coupled to the system by the so- 
called Brans-Dicke (BD) parameter to. However, the interest in BD theory in recent years 
was mainly motivated by the fact that the scalar-tensor part of the low energy effective 
action of superstring theory resembles BD theory [p||. 

Global topological defects in Brans-Dicke like theories have been studied extensively |L9|] . 
In spherically symmetric solutions of SU(2) Brans-Dicke- Yang-Mills-Higgs (BDYMH) 
theory in Schwarzschild-like coordinates have been constructed. This includes both the 
globally regular monopole as well as the corresponding black hole solutions. 

In this paper, we construct axially symmetric monopole solutions of SU(2) BDYMH 
theory. We give the Lagrangian, the Ansatz and the boundary conditions in Section II and 
present our numerical results in Section III. The conclusions are summarized in Section IV. 
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II. SU(2) BRANS-DICKE-YANG-MILLS-HIGGS (BDYMH) THEORY 



In [^0| it was argued that for conformally invariant fields such as the Yang-Mills fields, 
the Lagrangian in the Einstein frame equals that in the BD frame. In the following, all 
solutions - unless otherwise stated - are those obtained in the Einstein frame. 



A. The Lagrangian 



The Lagrangian C of Brans-Dicke theory reads [16] : 

i6ttG v m 



£ = 77-7; ( ^ - TA^* + £m (1) 



G denotes Newton's constant and uj > —3/2 is the Brans-Dicke parameter. For \u\ — > 00 
standard Einstein gravity is recovered. A conformal transformation to the Einstein frame 
yields the Lagrangian C : 

with \l/ given by the BD scalar field \l/ : 



The matter Lagrangian C m is given by : 

C m = —F^F^ - ^ e -^ V ^^D^ a D^ a - e -2~f^G*v(3> a ) , (4) 



with Higgs potential 



V($ a ) = ^($ a $ a -r/ 2 ) 2 , (5) 



the non-abelian field strength tensor 

F a , v = dpAZ - d v Al + ee ahc A\Al , (6) 

and the covariant derivative of the Higgs field in the adjoint representation (a = 1, 2, 3) 

D^ a = d^ a + ee abc A b ^ c . (7) 

A and 77 are the Higgs field's self-coupling constant and vacuum expectation value, respec- 
tively and e denotes the gauge coupling constant. The prefactors of the covariant derivative 
and the Higgs potential in (Q) result from the conformal transformation from the BD to the 
Einstein frame : 

Note that the Lagrangian (0) resembles that of Einstein- Yang-Mills-Higgs-dilaton theory 



22| , |23|| with a specific coupling of the scalar field. 
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B. Axially symmetric Ansatz 



The axially symmetric Ansatz for the metric in isotropic coordinates reads |24 
ds 2 = -fdt 2 + j (dr 2 + r 2 d6 2 ) + - f r 2 sin 
For the gauge fields we choose the purely magnetic Ansatz [[K] 



A t a = , A r a = —v, a 



er 



Aa a 



H- 



if 5 



Tt 

A/ = -- sin 9 (H 3 v r a + (1 - H±)v e a ) 



while for the Higgs field, the Ansatz reads ||10|| : 

$ a = n{^ x v r a + § 2 ve a ) ■ 

The vectors v r , vg and are given by: 

(sin 9 cos nip, sin 9 sin ny?, cos 9) 
(cos # cos rup, cos ^ sin ny, — sin i 
(— sin rup, cos n<p, 0) . 



(9) 



(10) 



(11) 



(12) 



(13) 



Here the winding number n enters the Ansatz for the fields. Since we are considering static 
axially symmetric solutions, the function /, /, m, Hi, H 2 , H 3 , if 4 , $i, $ 2 and \l/ depend only 
on r and 9. The spherically symmetric Ansatz for the construction of n — 1 solutions in 
isotropic coordinates is recovered if the dependence on 9 is dropped and additionally m = I, 
H t = H 3 = $ 2 = 0, H 2 = H 4 . 

We introduce the following dimensionless variable x and the dimensionless field ip : 



x = er/r , ip = V8nG^ 
With this rescaling, the set of differential equations arising from the variation 

$S = 5^J Cy/^d 4 x^j = 

depends only on u and the following dimensionless coupling constants : 



a: 



47T- 



eM. 



pi 



V2M ] 



(14) 



(15) 



(16) 



ir 



where M w = er\ is the gauge boson mass, M H = V2Xr] is the Higgs boson mass and 
Mpi = G^ 1 is the Planck mass. 

The energy E{n) of the solutions for this choice of Ansatz is given in terms of the 
derivative of the metric function / at infinity : 



E(n) 



2a 2 



lim x 2 d x f . 



(17) 
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C. Boundary conditions 



We have to impose 10 conditions on each of the boundaries x = 0, x = oo, 9 = and 
9 = § [^5| to solve the set of 10 second order partial differential equations. Regularity at 



the origin requires : 

d x f(o,e) = d x i{o,e) = d x m(Q,0) = o, 3^(0,0) = (18) 

Hi(0, 6) = 0, i = 1, 3, #(0, 0) = 1, i = 2, 4, $,(0, 0) = 0, i = 1, 2 (19) 

At infinity, the requirement for finite energy and asymptotically flat solutions leads to the 
boundary conditions: 

/(oo, 0) = l(oo, 9) = m{oo, 9) = 1, ^(oo,0) = O (20) 

Hi(oo,9) = : i = 1,2, 3, 4, $ 1 (oo,0) = l, $ 2 (oo,0) = O (21) 

In order to obtain the right symmetry for the solutions, we set on the z-axis as well as on 
the p-axis (9 = and 9 = ir/2, respectively) : 

H x = H 3 = $ 2 = (22) 

and 

d e f = d e m = d e l = d e H 2 = d e H 4 = d B $i = d e ^ = (23) 



III. NUMERICAL RESULTS 



We have solved the set of partial differential equations numerically. First, we studied the 
behaviour of the solutions for fixed u, j3 and increasing gravitational coupling a. In FIG. 1, 
we show the energy per winding number E/n for (3 = 0, n = 2 and three different values 
of a;, uj = oo corresponds to the pure Einstein gravity limit, studied in |L5| and for a = 0, 
the BPS multimonopoles are recovered with energy per winding number (in our rescaled 
variables) equal to unity. For uj < oo and a = 0, the energy per winding number is that 
of the BPS multimonopoles but rescaled by the prefactor (2u + 3)/ (2u + 4) of the matter 
Lagrangian C m (see (Q)). The reason for this is that for a = 0, the Brans-Dicke function 
and the metric functions are trivial tp = and / = m = I = 1, respectively. Thus the 
matter Lagrangian £ m , which is proportional to the energy density, is that of the flat space 
Yang-Mills-Higgs multimonopoles. This reasoning is, of course, also true for (3^0. Thus : 

2Ct^ I ^ 

E(a = 0, n, u, (3) = — — -E(a = 0, n, u = oo, /3) . (24) 
2oj -\- 4 

For a —>■ a max the branch of Brans-Dicke multimonopoles bifurcates with the branch of ex- 
tremal Reissner-Nordstrom (RN) solutions. The extremal RN solution in the model studied 
here is given by : 
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X 




, a = «\/^— TT ( 25 ) 
V x + an ) • 

and 

#i(z) = $ 2 (x) = , i = 1, 2, 3, 4 , $x(x) = l, tj}{x) = (26) 

Note that the horizon of the extremal RN solution in isotropic coordinates is located at 
x = Xh = 0. The functions of the limiting solution thus correspond to those of the RN 
solution on the full interval x e [0 : oof. The solution ( p5|) and (p6|) has magnetic charge P 
and energy i?(n) depending on the Brans-Dicke parameter uj : 



P = n l£rl' E(n) = P/a* (27) 

In FIG. 1, we demonstrate that for all chosen values of uj, the branch of globally regular 
monopole solutions bifurcates with the corresponding branch of extremal RN solutions. 
Moreover, we observe that a max (uj), the maximal value of a for which globally regular 
monopole solutions exist, is increasing drastically for decreasing uj (see also FIG. 3). We 
find, e.g. for (3 = 0, n = 2 : 

cw(oo) « 1.49 , cw(0) « 1.89 , a max {-l) « 2.80 (28) 



This agrees with the results in , where it was found that for the n — 1 solutions the ratio 
a m ax(0) / a m ax(oo) w 1.3. 

In FIG. 2, we demonstrate the behaviour of the metric function / in the limit a —>■ a max 
for uj = 0, n = 2 and (3 — 0. / is shown as function of the compactified coordinate 



z = x/(l + x) [26|. For increasing a, the value of the metric function at the origin, /(0), 
decreases to the RN value /i?jv(0) = 0. Moreover, the angle dependence diminishes indicating 
that the limiting solution is the spherically symmetric RN solution. For a — > a max ~ 1.89, 
/ approaches the function given by (p5|). 

For u = oo and a = 0, like-charged monopoles are either non-interacting (in the limit of 
vanishing Higgs boson mass) or repelling. When gravity comes into play, it is possible for 
the multimonopoles to reside in an attractive phase ||15|| . In the BPS limit (/3 = 0), this is 
not suprising, for (3 ^ however it is apparent that for small values of f3 < f3(n) gravity is 
able to overcome the repulsion of the long-range magnetic field. In [[TJ| it was found that 
(3{n — 2) « 0.21. Nevertheless, the extension of the attractive phase is limited by the fact 
that the solutions exist only up to an n- and (3- dependent maximal value of the gravitational 
coupling a. We find that like Einstein gravity Brans-Dicke gravity is able to overcome the 
repulsion between like-charged monopoles for sufficiently high values of the gravitational 
coupling a. In FIG. 3, we show a eq (/3,u), the value of a for which the mass of the n = 1 
monopole and the mass per winding number of the n = 2 multimonopole equal one another, 
as function of uj for two different values of (3. For both values of /3 — 0.1 and (3 = 0.2, a eq is 
increasing for decreasing ui from its value at uj = oo. Equally, the value a max is increasing 
with decreasing uj. The multimonopoles reside in an attractive phase in the parameter space 
above the a eg -curve and below the corresponding a maa; -curve. Thus in analogy to the uj = oo 
limit Brans-Dicke gravity is able to overcome the repulsion of the long-range U(l) field for 
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(3 = 0.1 and (3 = 0.2 and sufficiently high values of a. Since a max is increasing drastically for 
small u, the extention of the attractive phase for small uj is much bigger than in the uj = oo 
limit. However, we find that the order of magnitude of the value A = E(n = 1) — E(n)/n, 
which is an indicator for the strength of attraction between like-charged monopoles doesn't 
change significantly within the whole domain of attraction in the a-w-plane. This can 
be seen from the table below where we show A (a;) = E(n = 1, uj) — E(n = 2,u)/2 for 



(n 



2,u)/k, k = 2,3,4 



Table 1 



k 


2 


3 


4 


A(uj = oo) 


0.0051 


0.0024 


0.0022 


A(uj = -1) 


0.0052 


0.0033 


0.0018 



In FIG. 4, we show the difference between the energy of the n = 1 monopole and the n 
multimonopoles, A = E{n = 1) — E(n)/n, as function of uj for a = 1, (3 = and n = 2, 3, 4, 
respectively. 

A stays nearly constant for a large range of uj and decreases for decreasing uj. For uj — > 
uj = —3/2, A decreases to zero which is due to the fact that the prefactor (2uj +3)/ (2uj +4) is 
equal to zero for uj = —3/2 and thus the mass of the solutions itself vanishes. It can be clearly 
deduced from FIG. 4 that the strength of attraction between like-charged monopolos is 
smaller in Brans-Dicke theory than in pure tensor gravity theory. This can be compared with 
the results in [^7] and for two different Einstein- Yang-Mills-Higgs-dilaton (EYMHD) 
models. EYMHD theory - like Brans-Dicke theory - constitutes a theory of gravity in which 
the metric tensor has a scalar companion, in the case of EYMHD, the dilaton. While in E7 



the standard coupling of the dilaton in 4 dimensions was studied, the model studied in f28 
arose from dimensional reduction of an Einstein- Yang-Mills system in (4 + 1) dimensions 
| 29| . In both models the dependence of the strength of attraction on the dilaton coupling 
was studied. It was found that the value A was increasing for increasing dilaton coupling 
with a and the dilaton coupling not too close to their maximal values. Since in the model 
studied here, \/87tG'j(uj) (see (|3])) can be interpreted as the analog of a dilaton coupling 
and is increasing for decreasing uj and fixed G (i.e. fixed a), the strength of attraction is 
decreasing for increasing "dilaton" coupling. 

Apparently, the strength of attraction is increasing for increasing n which suggests that 
large clumps of monopoles might be possible. However, FIG. 4 also gives a hint that 
monopoles in such clumps will be more bound in pure Einstein gravity (uj = oo) than 
in Brans-Dicke gravity for equal gravitational coupling a. 



IV. CONCLUSIONS AND SUMMARY 

We have studied axially symmetric multimonopoles with n > 1 in Brans-Dicke gravity. 
We find that for a maximal value of the gravitational coupling a these solutions collapse to 
form an abelian black hole which outside the horizon is described by an extremal Reissner- 
Nordstrom (RN) solution with trivial Brans-Dicke scalar field. The RN solution has mass 
and magnetic charge depending on the Brans-Dicke parameter uj. This behaviour can be 
compared to that in Einstein- Yang-Mills-Higgs-dilaton (EYMHD) theory |^3|,^7|. When 
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either the gravitational coupling or the dilaton coupling reaches its maximal value, the 
branch of monopole solutions bifurcates with the branch of extremal Einstein-Maxwell- 
dilaton (EMD) solutions. These have a naked singularity and a non-trivial scalar dilaton 
field. 

In analogy to the n — 1 solutions, the maximal value of the gravitational coupling up to 
where the globally regular monopole solutions exist, is increasing for decreasing Brans-Dicke 
parameter. Equally, a eq , the value of a where the mass of the n = 1 monopole and the 
energy per winding number of the n = 2 multimonopole equal one another, is increasing 
for decreasing uj. At the same time, a max , the maximal value of the gravitational coupling 
a up to where globally regular multimonopole solutions exist, is increasing. Remarkable is 
that the rate of increase for decreasing uj is much bigger for a max than for a eq . Thus, the 
extension of the attractive phase is much bigger for small uj than for uj = oo. We conclude 
that at comparable values of the gravitational coupling a, the Brans-Dicke monopoles are 
less bound than the pure Einstein gravity monopoles, and though the former can exist to 
much bigger values of the gravitational coupling than the latter, the strength of attraction 
is of the same order of magnitude for all values of a and uj for which bound multimonopoles 
exist. These phenomena can be explained by noticing that in the gravitational field equations 



always the combination a e ff := a^J(2uj + 3)/(2uj + 4) appears. This can be interpreted as 
an "effective" gravitational constant. Clearly, for a fixed and uj decreasing from infinity, this 
expresssion is decreasing. This explains why the value a max up to where the Brans-Dicke 
multimonopoles exist is increasing for uj decreasing and also why at fixed a, the strength of 
attraction is decreasing for decreasing uj. Only at comparable values of a for different uj (e. 
g. at the k-th part of a max (uj)) is the strength of attraction of the same order of magnitude 
in Einstein theory and Brans-Dicke theory, respectively. 

In p0| the corresponding spherically symmetric black hole solutions were studied with 
emphasis on their thermodynamic properties. Axially symmetric black hole solutions in 
Einstein- Yang-Mills-Higgs theory have been studied recently |3(J in the context of the so- 



called "Isolated horizon framework" pTJ. It would be interesting to analyse in which sense 



the corresponding black hole solutions of the model studied in this paper fulfill the predic- 
tions of this framework. 
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FIG. 1. The energy per winding number E/n of the axially symmetric n = 2 Brans-Dicke 
monopole (solid) is shown for [3 = 0, u = oo, and —1 as function of a. Also shown is the energy 
per winding number of the corresponding Reissner-Nordstrom (RN) solution (dotted). 
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FIG. 2. The metric function / is shown as function of the compactified coordinate z = x/(l+x) 
for n = 2, to = 0, (3 = and several values of a including a = a max « 1.89. 
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FIG. 3. The value <x e q, where the mass per winding number of the n — 2 multimonopole and the 
n = 1 monopole equal one another is shown for (3 = 0.1 and (3 = 0.2 as function of u> — u>o, = — §■ 
Also shown is a max for n = 2 and the same values of (3. Note that the values of a max increase 
drastically with decreasing u, e.g. a max {f3 = 0.1, a; = —1) w 2.7 and a max (/3 = 0.2, cj = —1) ~ 2.6. 
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FIG. 4. The difference A = E{n = 1) — E{n)/n between the energy per winding number of the 
n = 1 monopole and the multimonopoles with winding number n is shown as function of uj — luq, 
ooq = — f , for a = 1.0, (3 = and n = 2, 3,4, respectively. 
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